Delocalised oxygen as the origin of two-level defects in Josephson junctions 
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One of the key problems facing superconducting qubits and other Josephson junction devices is 
the decohering effects of bi- stable material defects. Although a variety of phenomenological models 
exist, the true microscopic origin of these defects remains elusive. For the first time we show that 
these defects may arise from delocalisation of the atomic position of the oxygen in the oxide forming 
the Josephson junction barrier. Using a microscopic model, we compute experimentally observable 
parameters for phase qubits. Such defects are charge neutral but have non-zero response to both 
applied electric field and strain. This may explain the observed long coherence time of two-level 
defects in the presence of charge noise, while still coupling to the junction electric field and substrate 
phonons. 



Decoherence is currently a major limitation for super- 
conducting qubits and Josephson junction based quan- 
tum devices in general. An important source of decoher- 
ence stems from environmental two-level systems [TJ [2] . 
Recent experiments have even probed these defects di- 
rectly and shown that they are stable, controllable and 
have relatively long decoherence times themselves [5H7]. 
Little is known about the true microscopic nature of these 
defects, although many phenomenological theories ex- 
ist [8HT2]. We take a novel approach to the problem: 
starting from atom positions and species, motivated by 
ab initio and molecular mechanics methods. Using this 
approach we compute experimentally observed param- 
eters such as resonant frequency, defect-qubit coupling 
and response to strain, and find excellent agreement with 
experiments. We show that the quantum property of de- 
localisation of the oxygen atomic position in aluminium 
oxide naturally results in a model for two-level defects 
without the need for additional impurities. 

The existence of bistable defects in glasses and amor- 
phous solids in general is well known [15] . Amorphous 
insulating barriers (either in the form of Josephson junc- 
tions (JJ) or simply a native oxide) form an integral part 
of superconducting circuits, so it comes as no surprise 
that two-level systems (TLSs) are often considered to be 
an important source of noise in these circuits [TJ [2j [5]. 
The recent development of controllable qubit circuits 
(charge, flux or phase) has provided the opportunity 
to study so-called 'strongly coupled defects' [3] [6j [JJ. 
These defects have comparable resonance frequencies to 
the qubit circuit and coupling strengths and decoherence 
times large enough to allow coherent oscillations between 
qubit and TLS. Probing individual defects has promoted 
their bistable nature from hypothesis to observable fact 
as well as providing clues to their microscopic origin. 

There exists an embarrassment of riches in terms of 
theoretical models for such defects. Various phenomeno- 
logical models exist, including charge dipoles [8 , Andreev 
bound states [9], magnetic dipoles [10], Kondo impu- 
rities [11] and TLS state dependence of the JJ trans- 



parency [12] . Although detailed fitting of experimental 
data can place limits on these models [13] . they all have 
enough scope within their free parameters to explain the 
observed behaviour - rendering them presently indistin- 
guishable. 

To make concrete predictions, a detailed microscopic 
model of these defects is required. In this paper we con- 
sider the origin of defects to be within the amorphous 
oxide layer itself [15] , rather than assuming defects stem- 
ming from surface states [16] or the accidental inclusion 
of an alien species. A pertinent example defect is the 
oxygen interstitial in crystalline silicon. For an O de- 
fect in c-Si, the harmonic approximation for atomic po- 
sitions cannot be applied due to the rotational symme- 
try of the defect as oxygen delocalises around the Si-Si 
bond axis [17] . This forms an anharmonic system with a 
quasi-degenerate ground state, even in a "perfect" crys- 
tal. As many different spatial configurations can exist 
in the AlO^ amorphous junction, it is our premise that 
positional anharmonicity arises within voids in this layer. 
This yields TLSs with unique properties based solely on 
atomic positions and rotation in relation to the exter- 
nal electric field. Starting from this ansatz, we compute 
parameters which have been measured directly in ex- 
periments on TLSs, including: TLS resonant frequency, 
qubit-TLS coupling and TLS energy/strain dependence. 

To begin our investigation, a J J was modelled us- 
ing molecular mechanics and Density Functional The- 
ory (DFT). This approach has previously been used to 
model defects and dopants in aluminium oxide [I8j [19] . 
A representation of our system after amorphisation of the 
AIO1.25 barrier layer is displayed in Fig. [IJl Experimen- 
tal O/Al ratios have been shown to be highly dependent 
on fabrication processes 20 , 21 . therefore a representa- 
tive stoichiometry of AIO1.25 was chosen for simplicity. 
Various values of the oxide density were calculated and 
3.2 g/cm 3 was found to minimise the total energy of the 
system (Fig. [1J), which agrees with experimental val- 
ues [22] , This density is 0.8 times that of a common 
crystalline form of AI2O3 (Corundum). Using these sim- 
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FIG. 1. a) Depicts a J J with two aluminium slabs surround- 
ing an amorphous AIO1.25 barrier (aluminium: gray, oxygen: 
red), b) The projected radial distribution function G(r) using 
oxygen as a reference, c) Shows the total energy per atom of 
this structure as a function of oxide density. Fluctuations in 
the fine structure of b and c are due to finite box restrictions of 
the model, d) An illustration of the 2D oxygen delocalisation 
model. Aluminium atoms in gray, with the delocalised oxygen 
atom probability density shown for an example ground state 
distribution. 



illations, we compute the projected radial distribution 
function G(r) (¥\g.^p) of the resulting atomic positions. 
Both the Corundum peak (~ 1.97 A) and a broad dis- 
tribution (> 3 A) corresponding to the amorphous MO x 
layer are visible. Details of the precise parameters used in 
the DFT calculations can be found in the supplementary 
material. 

The energy scale for JJ defects observed in experiments 
is < 40 /ieV |H El 13- This energy splitting, while large for 
qubit experiments, is very small when compared to typ- 
ical electronic structure calculations, i.e. the ground and 
first excited state of our defect form a quasi-degenerate 
ground state on the scale of crystal defect energies, which 
puts them below the precision limits of DFT. Using the 
G(r) data obtained from DFT as a starting point, we de- 



velop an effective single-body model. This allows higher 
precision calculations as a function of atom locations, 
using empirical potentials for the interactions between 
an oxygen and nearest neighbour aluminium atoms. We 
initially consider a cubic lattice of six aluminium atoms 
with an oxygen atom delocalised at its center as our pro- 
totype defect (Fig. displays a representation of the 
case for delocalisation in 2D). As the experimental results 
point to bistable defects, we assume that the observed 
behaviour does not rely on a spherically symmetric po- 
tential in all three spatial dimensions, which would lead 
to triple degeneracy. We therefore concentrate on two- 
dimensional delocalisation. Using the empirical Streitz- 
Mintmire potential [23] we derive an effective single par- 
ticle Hamiltonian 
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h 2 



2m n 
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where m oxy is the mass of an oxygen atom and V(r) 
is the potential due to the six aluminium atoms, gen- 
erated by the Streitz-Mintmire formalism. The result- 
ing single-body time-independent Schrodinger equation 
is then solved on a finite grid (see supplementary mate- 
rial). 

For our effective model, the atomic positions |X|, \Y\ 
and \Z\ labeled in Fig.[l]i represent aluminium atom pairs 
(e.g. — X,+X) lying on the cardinal axes and displaced 
equidistantly from the origin in each direction. 

In two dimensions, the potential landscape of interest 
approximates the venerable "Mexican hat" potential as 
this is the rotationally symmetric generalization of the 
double-well potential in ID. This potential has a unique, 
spherically symmetric ground state and a doubly degen- 
erate first excited state, which we see for \X\ = \Y\ > \Z\. 
Small deformations due to \X\ 7^ \Y\ or translations of 
\X\ or \Y\ off axis quickly result in a quasi-degenerate 
ground state of the form seen experimentally. The re- 
sulting splitting, £?oi , is plotted in Fig. |2^i for a range of 
\X\ and \Y\ values, with a fixed value of \Z\ = 2.5788 A. 
Fig. shows the absolute dipole moment response over 
the same phase space. Subfigures [2]>f display four posi- 
tions of equal energy on Fig. where the first excited- 
state wavefunction of the oxygen atom is plotted with its 
effective potential. 

The phase space in figures |2^l and [2]o are split into 
four domains, [2^i via dispersive peaks with their max- 
ima located at the bifurcation lines visible on (2)3. The 
properties of these domains can be explained through the 
interplay of potential configuration and dipole alignment. 
In two-dimensions these potentials can be described as a 
sum of two ID potentials: one in the x direction and 
the other along y. The first case is a sum of two dou- 
ble wells (tetra-well) and the second, a sum of a double 
and harmonic well (hemi-tetra-well). Points c and d in 
Fig. [2] are examples of the former while e and f represent 
the latter. Similarly, each domain has a dipole element 



3 




»1H 

-4.5 -2.25 2.25 4.5 -4.5 -2.25 2.25 4.5 

x (A) x (A) 

b'b 



X(k) 



X (A) 



FIG. 2. a) Map of the Eoi energy splittings of the delocalised oxygen 2D model. The \X\ and \Y\ axes represent aluminium 
pair positions with \Z\ = 2.5788 A. b) The difference between the absolute dipole moment (in x- and ^/-directions) over the 
same range. We see either \p x \ (red) or \p y \ (blue) dominated behaviour in all regions except \X\, \Y\ < 3.5 A (where the 
oxygen is tightly confined), c) through f) show the first excited state wavefunction of the oxygen atom and the acting potential 
of four configurations indicated on a). For comparison with existing qubit experiments we plot contour lines corresponding 
Eoi/h = 0.5 — 10 GHz (red to yellow) overlayed on a) and b). The same resonant frequencies are discussed in Fig. [3] hence c-f 
all represent configurations of Eoi/h = 8 GHz; whereas points i-vi are the locations selected for the strain study in Fig. [4] 



which is orientated in either x or y, as one observes from 
the direction of the nodal line in each subfigure. These 
computed dipole moments correspond well to observed 
values, assuming O (nm) junction widths [8| fl4] . 

To compare our TLS model directly to experiments, we 
assume that our J J lies within a phase qubit, although 
the model applies equally for any device comprised of 
amorphous junctions. The measurable signal of a TLS 
in a phase qubit is the resonance of the TLS and qubit 
splitting energy, £?oi, with the qubit-TLS coupling, S max . 
For the phase qubit [8], S max is a function of Eoi and 
p [24], the effective dipole moment due to an electric 
field applied in the direction of delocalisation, 
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E, 
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(see supplementary material). Throughout this discus- 
sion we assume a junction width w = 2 nm and capaci- 
tance C = 850 £F. 

In Fig. [3] we plot contour lines representing constant 
values of £"01 which correspond to the purview of exper- 
imentally observed qubit resonant frequencies. This re- 
gion of parameter space corresponds well with the calcu- 
lated G(r) (see Fig. [Ip) in that spacings of \X\ = 4 — 6 A 
are possible, albeit uncommon. This graph shows the 
\X\ response in the domains where \p x \ is dominant, al- 
though it is clear to see that the |Y| domains are sym- 
metric from Fig. [2] The S max (Eq. |2| response to these 
frequencies is plotted as a function of \X\, in which we 
see maximum coupling strengths which correspond well 
with experimental observations [5j El [14] . Two contours 
are visible for each £qi splitting, which is due to the bi- 
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FIG. 3. Coupling strength to a fictitious phase-qubit Smax 
as a function of \X\ in the domains where \p x \ is dominant 
(see Eq. [2| for a set of constant £01 splitting frequencies. 
For comparison with experimental results, £01 and Smax are 
expressed in frequency units. 



furcation in Fig. [2] as the values of \X\ and \Y\ flip in 
phase space - indicating a change from the tetra- to the 
hemi-tetra region. 

A key observation of the TLS-qubit experiments is the 
unusually long coherence times of strongly coupled de- 
fects [4] [25]. As our model assumes a charge- neutral 
defect, coherence time is linked to the dipole element 
(for charge noise) and the strain response (for phonons). 
The strain response has recently been observed directly 
through mechanical deformation of a phase-qubit [26 . 

We introduce a series of deformations in our 2D model 
to measure the variation in £01, which are depicted in 
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FIG. 4. a) depicts a number of deformations which were ap- 
plied to the aluminium atoms in the x — y plane (see text). 
Only one of the modes responded with over a few Hz of move- 
ment (highlighted), which is indicative of an optical phonon 
mode, generating frequency splittings of O (100 MHz) for pi- 
cometer deformations. A deformation range of 100 pm was 
applied to the points i-vi from Fig. [2] yielding responses that 
are both hyperbolic and symmetric (b). Over the range 
20 — 100 pm the response is linear with a strain gradient, 
shown in c) for both the tetra- and hemi-tetra domains in the 
^-direction. 



Fig. [4^i. All deformations were tested in each of the four 
regions of Fig. [2^i, not only in the x-direction as shown, 
but also in y. Of the tested deformations we find the 
response of one (the optical phonon mode, highlighted in 
Fig. [4^i) to be 10 5 times stronger than the others. Such 
a deformation corresponds to a translation of both alu- 
minium atoms in the same direction and relative to the 
oxygen, along the axis of the dominant dipole. This sug- 
gests an explanation for the long TLS coherence times, 
as a delocalised oxygen is only sensitive to a small subset 
of available phonon modes (as well as coupling to charge 
noise only through its electric dipole). Fig. shows 
this response for the optical mode (at several points of 
interest labeled in Fig. [2^), displaying a characteristic hy- 
perbolic response which is typical of a two-level system. 
This compares well with the observed strain response in 
Ref. [26| Finally, Fig. ^ shows the linear strain gradient 
plotted along the E^i/h = 8 GHz contour for the tetra- 
and hemi-tetra- regions in the \p x \ ^ 0, \ p y \ =0 domains. 

Our model allows prediction of experimentally mea- 
sured properties of strongly coupled TLSs with atomic 
positions as the only input parameters. Using realistic 
atomic positions obtained via molecular mechanics and 
ab initio methods, the correspondence with observed de- 
fect properties is excellent and therefore suggests that 
these defects can arise in AlO^ without any alien species 



present. Our model also proposes that restricting the de- 
localisation of oxygen, for example through higher den- 
sities in the amorphous layer, results in fewer voids [15] 
and therefore fewer two- level defects. Microscopic mod- 
els of this type will guide future fabrication and design of 
superconducting circuits, leading to lower levels of noise 
and greater control over their quantum properties. 
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